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SECTION  1 


INTRODUCTION 


This  report  is  intended  to  address  three  questions  that  arise  during  simu¬ 
lation  or  hardware  testing  activities  of  communications  links  under  Rayleigh  fading 
conditions:  How  many  decorrelation  times  per  realization  are  necessary?  How  many 
samples  per  d  correlation  time  are  necessary?  How  should  interpolation  be  done  be¬ 
tween  samples?  These  questions  are  answered  in  part  in  the  DMA  signal  specification 
for  nuclear  scintillation  (Wittwcr  1980)  which  requires  a  minimum  of  100  decorrela¬ 
tion  times  per  realization  and  10  samples  per  decorrelation  time.  However,  experience 
has  shown  that  receiver  performance  will  show  considerable  statistical  variation  when 
the  minimum  realization  length  is  used.  This  is  particularly  true  of  links  which  have 
large  power  margins  and  arc  susceptible  to  only  t  he  very  deepest  fades.  Of  course  the 
best  way  to  answer  these  questions  is  to  measure  link  performance  with  realizations  ol 
increasing  length  and  resolution  until  the  statistical  variation  in  the  results  from  one 
realization  to  the  next  is  acceptable.  Unfortunately,  the  luxury  of  doing  this  analysis 
often  docs  not  exist. 

The  next  higher  level  of  analysis  of  these  questions  is  to  look  at  the  statistics 
of  the  realizations.  This  is  the  approach  that  will  be  taken  in  this  report.  The 
first  order  statistics  of  the  realizations  arc  measured  by  calcu'ating  the  cumulative 
distribution  of  the  amplitudes  and  comparing  this  with  the  Rayleigh  distribution. 
The  second  order  statistics  of  the  realizations  are  measured  by  calculating  the  mean 
duration  and  separation  of  fades  and  comparing  these  quantities  to  ensemble  values 
for  Rayleigh  fading. 

In  general,  the  received  signal  may  be  written  as  the  convolution  of  the 
channel  impulse  response  function  h(t,r)  with  the  transmitted  signal  s (/.): 


roc 

u(t)  -  /  h(t,T)s{t  -  t  )dr  (LI) 

J  o 

In  either  software  link  simulations  or  in  hardware  channel  simulators,  Kquat.ion  1.1 
can  be  implemented  as  a  tapped  delay  line: 

A/>  J 

u(0  77  /f(',yAr).s(<  -jAt)At  (1.2) 

)  o 
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wlim.  N n  is  number  ('['  taps  on  tlic  delay  line;  A r  is  the  delay  spacing  of  the  delay  line; 
h(t,j Ar)  is  the  time  varying  complex  weight  of  the  jl/*  tap;  and  s(t)  is  the  input  signal. 
In  a  software  simulation  of  link  performance,  time  will  also  be  discretely  sampled  (i.e. 
I  --  kAt). 

Under  Rayleigh  fading  conditions,  /i(t,r)  is  a  complex,  zero  mean,  normally 
distributed  random  variable  and  thus  has  a  Rayleigh  amplitude  distribution.  It  then 
follows  from  Equation  1.2  that  u(<)  is  also  a  complex,  zero  mean,  normally  distributed 
random  variable  with  a  Rayleigh  amplitude  distribution. 

A  complete  analysis  of  these  issues  would  consider  the  sampling  requirements 
for  each  delay  of  the  discrete  impulse  response  function  /t(ArAf,  j  Ar).  However,  this  is 
beyond  the  scope  of  this  report.  Therefore,  satnplmg  requirt  ments  on  the  flat  fading 
impulse  response  function  h[kAt ),  where 

;V„.  1 

h(kAl)  ^>2  h[kAt,jAr)AT  ,  ( 1 

;  o 

will  be  addressed  in  this  report.  The  sampling  requirements  for  h(lcAi)  will  give  some 
indication  of  the  sampling  requirements  for  the  frequency  selective  impulst  response 
function  h(kAt,j  At).  Perhaps  this  should  he  stated  another  way:  Sampling  that  is 
inadequate  for  h(kAt)  wili  surely  be  inadequate  for  h(kAt,j  Ar).  Thus  it  is  the  intent 
of  this  report  to  define  adequate  sampling  for  h(kAt)  and  to  inter  adequate  sampling 
requirements  for  each  delay  of  h(kAi,j  At). 
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SECTION  2 


temporal  statistics  of  scintillation 


This  section  is  a  review  of  well  known  results  from  the  classical  work  ol 
S.  ().  Rice  (1918,  19.T1,  I9a8)  on  the  first  and  second  order  statistics  of  Rayleigh 
fading. 


2.1  Fill  ST  ORDER  STATISTICS. 


Under  strong  scattering  conditions,  the  electric  held  incident  on  the  plane 
of  the  rereiver  is  the  summation  of  many  waves  propagating  in  slightly  different  di¬ 
rections  a  Lout  the  imo-of-sight.  Under  the  centre,!  limit  theorem  of  statistics,  the 
two  orthogonal  components  of  the  electric  field  must  then  he  zero-meati,  normally 
distributed  random  variables.  It  is  assumed  that  the  two  orthogonal  components  are 
aLo  independent.  The  complex  narrow-band  envelope  of  the  electric  held  undergoing 
Rayleigh  fading  may  be  then  represented  as 


K{()  ■  *(0  !  «J/(0 


(2.1) 


where  i  and  y  are  independent  and  normally  distributed  with  v.ero  mean  standard 
deviation  o.  The  carrier  frequency  term,  exp(tuV),  has  been  neglec  ted  in  this  expres¬ 
sion.  Thus  /,’(/)  may  be  thought,  of  as  the  output,  voltage  of  a  down-converter  where 
x(t)  is  the  in-phase  component  and  y(l)  is  the  quadrature-phase  component. 

It.  is  easy  to  show  that  the  probability  density  fuuciioii  of  the  amplitude  ol 


is  Rayleigh: 


«(0  0  y-(') 


(2.2) 


f[“)  ,,  'M> 

‘  t* 


(2.: 
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where  PQ  is  the  mean  power  (P0  -  (a2)  -  2 o2  where  the  brackets  (•}  denote  on 
ensemble  average)  The  cumulative  distribution  of  the  power  P ,  which  is  equal  to  the 
probability  that  the  power  is  less  than  oi  equal  to  P,  is  given  by 


rVr 

P 

/  J  (a)da  —  1  —  exp 

Jo 

■  Pa 

(2-1) 


This  well  known  function  is  plotted  in  Figure  2.1  versus  the  ratio  P/ P0.  It  can  be  seen 
from  the  figure  that  the  piobability  of  a  20  dB  fade  or  deeper  below  the  mean  power 
level  is  10~2  and  the  probability  of  a  30  dB  fade  or  deeper  is  10-3.  This  is  consequence 
of  the  fact  that  for  smali  values  of  P/Pa ,  F{P)  ~  P  /  P0. 


Figure  2.1.  Cumulative  distribution  of  Rayleigh  fading. 
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2.2  SECOND  ORDER  STATISTICS. 


The  fading  rate  is  determined  by  the  second  order  statistics  of  the  electric 
field.  The  autocorrelation  of  the  electri .  fief1  is,  in  general, 


{K(t)E'(t  -r  r))  =  {x(t)x{t  F t))  -i  (. y{t)y[t  +  r))  -  2o2p(t)  .  (2.5J 

There  are  two  limiting  forms  for  the  autocorrelation  function  p(r).  Under  strongly  dis¬ 
turbed  scattering  conditions  th  *  occur  at  early  times  or  at  the  center  of  the  disturbed 
region,  p(r)  has  the  Gaussian  form 


p{r)  -  exp 


(2 .6) 


where  Tq,  the  decorrelation  time  of  the  electric  field,  is  defined  as  the  e  folding  point 
of  the  autocorrelation  function  (p(r0)  -  e The  corresponding  Doppler  spectrum 
of  the  temporal  fluctuations  is 


j 


f  oo  . 

/  exp(  — Jcur )p(r )d.T  -  V7rr0exp 


(2.7) 


which  also  has  the  Gaussian  form.  Under  kiss  highly  disturbed,  but  -’till  Rajleig.i 
fading  conditions,  the  autocorrelation  function  is  usually  assumed  to  have  the  form 


1  1 - 

exp 

a  '■/' 

7o 

To 

(2.K) 


where  the  parameter  o(a  --  2.1iGlt)2)  is  determined  by  the  ( oi  hi.ion  ■  li.it  c 

'Jdie  corresponding  Doppler  spectrum  has  the  form  commonly  reierred  to  as  an 
spectrum: 


J  In 

-  .1  • 


(2.‘>) 


A  comparison  of  realization.,  of  the  impulse'  rcM’tm.-.f  funclion  w.tli  Gaus¬ 
sian  and  f  4  Doppler  speetrums  are  sliown  in  Figure  2  wi  t  the  signa  power  in  dll 


is  ploU'ul  versus  normalized  time/r0.  These  realizations  have  been  generated  using 
statistical  techniques  described  in  Appendix  H.  both  realizations  were  generated  from 
the  same  set  of  random  numbers  so  there  is  a  strong  correlation  in  the  features  seen 
in  the  two  realizations.  Also,  both  realizations  have  unity  mean  power.  The  f  4  re¬ 
alization  has  a  more  spiky  appearance  due  to  having  more  energy  at  higher  Doppler 
frequencies.  The  two  signals  have  similar  low  frequency  behavior  and  the  fades  of  the 
two  signals  follow  each  other  quite  closely.  The  difference  in  the  two  signals  is  the  high 
frequency  jitter  of  the  f  4  signal  about  the  more  smoothly  varying  Gaussian  signal. 
The  significance  of  this  on  the  temporal  statistics  of  the  fades  wjl!  become  apparent 
lat  er. 


2.3  TEMPORAL  STATISTICS  OF  RAYLEIOH  FA  DIN  O . 


The  mean  duration  and  separation  of  fades  below  an  arbitrary  nower  level 


r  aim  i  n a i  oi  nares  a on\  »• 
crossings  of  the  level  I’  in  the  time  interval  T. 
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The  probability  that  the  amplitude  a  will  cross  the  level  f.  d  l‘  in  the  time 
interval  t  to  t  •.  dt  with  a  positive  derivative  is  equal  to  the  probability  that  a  >  0 
and  that  f  a' dt  ■'  a  •;  P.  This  prohabiiity  is  given  by 


/(«,«*)  dt  (  da!  af(f,a) 
J  o 


where  /(»,«')  is  the  joint  probability  density  function  of  t  he  amplitude  a  and  its  time 
derivative  a1  -  dn/ilt.  The  probability  that  a  will  cross  the  level  P.  in  the  time  interval 
Mot  I  dt  with  a  derivative  of  either  sign  is  then 

dt  [  \u'J[P,<i)dd  . 

J  Ou 


lor  stationary  processes,  the  mean  number  of  level  crossings  of  I ’  in  the  interval  t  to 
t  ■  T  I  hen  becomes 


[S\l\T)i  T  r  I r,',/(/\«')V 

J  (X. 


(2.10) 


The  join  t  probabi  I  it  y  ( leu  si  l  \  fu  net  ion  of  the  am  pi  it  ude  ami  il  s  t  i  i  ne  deriva¬ 
tive  is  (Kiee  1018,  Dana  1082) 


2 -a 


a7  ,  i  r0/&  \ 

M'/A)2] 
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2  P0 
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(0  <  a  <  oo,  -oo  <  a'  <  oo) 

This  probability  density  function  is  derived  in  Appendix  A.  It  can  be  seen  from  the 
form  of  Equation  2.11  that  the  probability  density  function  of  a  is  Rayleigh;  that  the 
probability  density  function  of  a'  is  normal  with  zero  mean  and  with  A2/o/tq  variance; 
and  that  a  and  a'  are  independent.  The  functional  form  of  f(a,a')  is  the  same  for 
either  of  the  limiting  forms  of  the  Doppler  spectrum.  The  only  difference  is  the  value 
of  the  parameter  A  which  depends  on  the  functional  form  of  the  Doppler  spectrum 
(A  =  1  for  the  Gaussian  spectrum  and  A=  1.518  for  the  f~4  spectrum). 

The  mean  number  of  level  crossings  can  now  be  easily  evaluated  with  the 

roc  nit 

i  ^  V*  »  V 


(N(P,T))  =  A 


[W 

*Po 


exp 


P 

To. 


(2.12) 


The  effect  of  the  different  Doppler  spcctrums  is  to  scale  the  mean  number  of  level 
crossings  by  the  quantity  A.  Recall  that  this  fact  was  shown  qualitatively  by  com¬ 
paring  the  two  random  realizations  in  Figure  2.  By  noting  that  two  level  crossings 
are  required  to  define  the  beginning  and  ttie  end  of  a  fade  below  I\  the  mean  rate  of 
fades  below  the  level  P  is  (A '(P,  r0))/2r0. 

Figure  2.3  shows  plots  of  the  mean  number  of  crossings  of  P  in  one  decorre¬ 
lation  time  versus  the  ratio  P/P(,  for  the  two  limiting  forms  of  the  Doppler  spectrum. 
The  maximum  value  of  ( N(P,t0 ))  occurs  at  P j P0  =  1/2  or  -3  dB.  As  a  point  of  com¬ 
parison,  the  median  level  M  of  the  Rayleigh  distribution  is  M  —  Po?-n( 2)  or  1.0  dll 
below  the  mean  power. 

'Fhe  mean  separation  ('/’<,■  (/'))  of  lades  below  P  is  obtained  from  the  mean 
number  of  fades  per  unit  time  ft  =  (N (P,  r0))/'2r0.  For  any  long  time  interval  T ,  the 
mean  number  of  lades  is  ftT  and  fhe  mean  separation  is  just  f  j  [VP  or  i/ft.  Thus  the 
mean  separation  of  fades  below  P  is 


<Ts(P)) 


2  jn  r" 
A  V  2  T 


exp 


(2.13) 
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Figure  2.3.  Mean  number  of  level  crossings  per  r0. 

The  mean  separation  of  a  fade  below  P  is  equal  to  the  average  time  between  crossings 
of  P  with  either  a  negative  value  of  a'  (which  defines  the  start  of  the  fade)  or  with 
positive  value  of  a'  (which  defines  the  end  of  the  fade).  Thus  the  mean  separation  of 
fades  below  P  is  also  equal  to  the  mean  separation  of  flares  above  P. 

The  mean  duration  (TfAde{P))  of  fades  below  P  is  obtained  as  follows. 
During  a  long  time  interval  T,  the  total  time  that  the  power  will  be  b  'low  P  is  F[P)T 
where  F(P)  is  the  cumulative  distribution  given  in  Equation  2.4.  The  mean  duration 
is  then  the  sum  of  all  durations  F(P)T  divided  by  the  number  of  fad'1'7  0T.  The  result 
is 


{Tfade[P)) 

TO 


(2.14) 


The  mean  duration  of  a  flare  above  P  is  the  mean  time  that 

the  power  stays  above  P.  Using  the  arguments  given  above,  the  mean  separation  of 
a  fade  or  a  flare  is  equal  to  the  mean  time  that  the  signal  is  above  P  plus  the  mean 
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time  that  it  is  below  P  :  (TfADe[P))  +  { Tflare{P ))  =  (7s(P)).  The  mean  duration 
of  a  flare  is  then 


{Tflare{P))  _  1  l*P°  / o  i  r \ 

- 7 - "A VTF  ■  (2'15) 

The  mean  duration  and  separation  of  fades  and  flares  are  shown  in  Figure  2.4 
for  a  Gaussian  Doppler  spectrum.  For  the  f-4  Dopplei  spectrum,  all  three  curves  in 
Figure  2.4  scale  b)  1/A  =  0.6589. 


Figure  2.4.  Mean  duration  and  separation  of  fades  and  flares. 
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SECTION  3 


SAMPLED  RAYLEIGH  FADING 


The  requirements  on  the  sampling  of  Rayleigh  fading  are  given  in  the  DNA 
signal  specification  for  nuclear  scintillation  (Wittwer  1980)  which  requires  a  minimum 
of  100  decorrelation  times  per  realization  and  10  samples  per  decorrelation  time.  The 
questions  that  arise  from  this  requirement  can  be  summarized  as:  How  close  are 
such  realizations  to  Rayleigh  fading'  To  address  this  question,  random  realizations 
of  Rayleigh  fading  will  be  generated;  moments  of  the  amplitude,  cumulative  distri¬ 
bution,  and  mean  fade  duration  v/ill  be  measured;  and  these  measured  values  will 
be  compared  with  their  ensemble  values.  Of  course,  because  of  the  finite  number  of 
samples  in  each  realization,  each  of  these  measurements  is,  in  fact,  a  random  variable 
with  some  mean  and  standard  deviation.  The  variations  in  these  measurements  from 
realization  to  realization  are  measured  by  generating  a  large  number  (102 T  to  be  ex¬ 
act)  of  realizations  and  computing  the  mean  and  standard  deviation.  A  method  of 
generating  such  realizations  is  outlined  in  Appendix  B. 

The  objectives  of  this  section  are  to  present  the  means  and  standard  de¬ 
viations  of  such  measurements,  and  to  attempt  to  answer  the  above  question  based 
or:  these  results.  The  effects  of  interpolation  between  realization  sampies,  and  the 
effects  of  using  longer  realizations  (more  decorrelation  times  per  realization)  will  also 
be  examined. 


3.1  MEASURED  STATISTICS  OF  100  r0  REALIZATIONS. 

Realizations  generated  with  N  — 102-1  time  samples  and  ;V0  .  10  samples 

per  decorrelation  time  wiii  serve  ;us  a  baseline  calculation.  These  realizations  are  also 
sampled  at  r0/H)  rime  intervals. 

One  criterion  fur  deciding  that  a  realization  has  Rayleigh  amplitude  statis¬ 
tics  is  that  the  moments  of  the  amplitude  .d.ould  agree  with  Rayleigh  values.  Ensemble 
values  tor  the  moments  of  the  amplitude  are  easily  obtained  from  Equation  2. -5  and 
are: 
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(3.1) 


(a) 

=  \\!*r c 

(a2) 

=  Po 

(a3) 

=  l\/nPo 

(O 

=  2/tf  . 

The  .scintillation  index  54  is  the  standard  deviation  of  the  power: 


S4 


(a-)  ~  (a2)2 
(a2)2 


(3-2) 


It  is  necessary,  but  not  sufficient,  that  St  equal  unity  for  Rayleigh  fading.  The  scin¬ 
tillation  index  is  a  good  measure  of  the  statistics  of  flares  but  net  of  fades.  Another 
statistic  that  is  sensitive  to  the  distribution  of  fades  is  the  mean  log  amplitude  x: 


X  {in  a)  =  in\fF0  -  h/2  (3-;i) 

where  7  is  Euler’s  constant  (7  --=  0.5772157- •  •). 

The  measured  mean  and  standard  deviation  of  the  amplitude  moments,  S’4, 
and  x  are  presented  in  Table  3.1.  The  measured  values  of  a  single  realization  should 
equal  the  ensemble  value  plus  or  minus  <  e  or  two  standard  deviations.  It  can  be 
seen  from  the  table  that  the  a\eragc  values  of  the  moments  are  close  to  the  ensemble 
values  but  that  the  standard  deviations  of  the  higher  moments  and  X  can  be  as  large 
as  25  percent. 

Perhaps  a  better  criterion  for  the  first  order  statistics  is  close  agreement 
between  the  Rayleigh  cumulative  distribution  and  the  measured  cumulative  distri¬ 
bution.  The  measured  cumulative  distribution  (dots  plus  or  minus  one-sigma  error 
bars)  is  plotted  in  Figure  3.1  along  with  the  ensemble  curve  (Kqn.  2 .4).  A  level  c.  0 
d I i  corresponds  to  the  mean  power  P0.  It  can  been  seen  from  the  figure  that  100  r0 
realizations  do  indeed  have,  on  trie  average,  a  Rayleigh  distribution  of  fades  down  to 
at  least  30  dB.  It  is  clear,  however,  that  the  possiole  deviation  from  Rayleigh  fading 
of  a  single  realization  bec<  s  larger  as  mie  examines  deeper  fades. 
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Tabic  3.1.  Amplitude  moments  of  100  r0  realizations. 


Moment 

Ensemble  Value 

Average*  Standard  Deviation* 

(a) 

y/fft/2 

1.003 

- - - — - 

0.057 

(a2) 

Po 

1.005 

0.112 

(a3) 

3\A^o/4 

1.007 

0.172 

(a4) 

i 

2/J02 

1.009 

0.243 

1.0 

0.983 

0.U83 

X 

-1/2 

0.991 

_ 

0.221 

*  Normalized  to  the  ensemble  values. 

The  fidelity  of  the  realizations  in  reproducing  the  second  order  statistics  of 
the  fading  will  be  investigated  by  considering  the  mean  fade  duration  and  separation. 
The  mean  fade  duration  is  a  good  statistic  to  examine  for  communications  applications 
because  errors  often  occur  in  bursts  during  deep  fades.  If  the  fades,  on  the  average, 
arc  too  long  or  too  short,  the  error  bursts  will  not  have  the  proper  durations  and  the 
resulting  receiver  performance  may  be  misleading.  Fade  duration  measurements  and 
the  ensemble  curve  for  a  Gaussian  Doppler  spectrum  are  shown  in  Figure  3.2.  because 
the  realizations  are  generated  with  10  samples  per  decorrelation  time,  the  minimum 
fade  duration  is  ro/10.  The  measurements  reproduce  ensemble  values  for  fades  down 
to  about  20  dll.  Below  this  level,  the  average  of  the  measurements  is  limited  to  a  value 
of  about  hi/ 10.  The  measured  value  for  30  d li  fades  is  less  than  r0/10  because  some 
100  r0  realizations  do  not  have  30  dll  fades.  The  effect  of  these  realizations  on  the 
measurements  is  to  reduce  the  average  values.  The  measured  and  ensemble  mean  fade 
separation  are  shown  in  Figure  3.3.  It  is  likely  that  100  r0  realizations  will  not  have 
two  fades  below  the  25  or  30  dll  level  so  a  fade  separation  measurement  can  not  he 


Figure  3.3.  Mean  fade  separation  of  100  t0  realizations. 

made.  Thus  the  error  bars  at  these  levels  are  large  and  the  measured  separation  for 
30  dB  fades  is  low  by  a  factor  of  almost  10. 

For  fades  less  than  20  dB  or  so,  the  100  t0  realizations  reproduce  the  ensem¬ 
ble  values  for  fade  probability,  tade  duration  and  fade  separation.  If  deeper  fades  are 
of  concern,  however,  then  100  r0  realizations  sampled  at  t0/10  may  not  be  adequate 
for  producing  accurate  error  rate  predictions.  Two  possible  alternatives  are  discussed 
below:  interpolation  and  longer  realizations. 


3.2  EFFECTS  OF  INTERPOLATION  AND  LONGER 
REALIZATIONS. 


Although  the  realizations  of  the  impulse  response  function  are  usually  gen¬ 
erated  with  samples  spaced  at  r0 / 10,  the  sampling  period  of  the  simulated  or  tested 
receiver  can  be  much  smaller  than  r0/lO.  There  are  at  least  three  approaches  to  this 
problem.  One  apptoach  is  to  sample  the  'unpulse  response  function  at  a  rate  equal  to 


the  sample  rate  of  the  receiver.  However,  if  the  sample  period  is  much  smaller  than  r0, 
this  results  in  very  large  realizations  if  each  realization  must  have  100  decorrelation 
times.  Another  approach  is  to  hold  the  impulse  response  function  constant  during 
the  period  r0/10,  and  change  it  abruptly  at  the  end  of  the  period  to  the  next  value 
of  the  impulse  response  function.  This  is,  in  fact,  the  approach  used  to  measure  the 
temporal  statistics  of  the  realizations,  and  in  principle,  this  procedure  is  acceptable 
because  t0/  10  sampling  should  result  in  small  changes  from  one  value  of  the  impulse 
response  function  to  the  next.  A  third  approach  is  to  interpolate  between  values  ot 
the  impulse  response  function.  Simple  linear  interpolation  will  be  considered  here. 

Figure  3.4  shows  the  mean  fade  duration  measured  from  exactly  the  same 
realizations  used  to  generate  the  results  of  the  previous  section.  However,  here  the 
realizations  are  sampled  with  a  period  of  ro/40  using  linear  interpolation  between  the 
ro/ 10  values.  During  the  r0/40  period,  the  impulse  response  function  is  held  constant. 
The  feature  to  note  in  Figure  3.4  is  that  the  measured  mean  fade  durations  now  lie 
on  top  of  the  theoretical  curve  for  fades  down  to  30  dB.  Of  course,  if  the  mean  fade 


Figure  3.4.  Mean  fade  duration  of  100  r0  realizations  with  interpolation. 
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duration  of  35  dB  fades  had  been  measured  it  would  be  equal  to  about  r0/dO  and  would 
lie  above  the  theoretical  curve.  Because  the  variations  in  the  mean  fade  duration  from 
realization  to  realization  are  determined  by  the  number  of  decorrelation  times  in  the 
realizations  rather  than  the  sampling  period,  the  variations  around  the  mean  values 
are  about  the  same  size  as  measured  without  interpolation. 

Finally,  the  cumulative  distribution,  mean  fade  duration,  and  mean  fade 
separation  of  400  tq  realizations  generated  with  10  samples  per  7o  (N  =  d096,  iV0 
—  10)  and  sampled  at  ro/dO  using  linear  interpolation  are  shown  in  Figures  3.5,  3.0, 
and  3.7.  Here  the  variations  around  the  mean  values  are  smail  compared  to  variations 
of  100  r0  realizations.  Also,  the  mean  separation  of  30  dl3  fades  for  dOO  7n  realizations 
is  close  to  the  theoretical  value  indicating  that  most  dOO  r0  realizations  have  at  least 
two  fades  of  this  depth. 

3.3  SUMMARY. 

It  is  clear  that  100  r0  realizations  do  not  accurately  reproduce  the  duration  or 
separation  or  fades  below  about  20  dB.  Linear  interpolation  can  be  used  to  accurately 
reproduce  the  duration  of  20  to  30  dB  fades  in  100  r0  realizations  but  it  takes  dOO  r0 
realizations  to  accurately  reproduce  the  mean  separation  of  30  dB  fades.  It.  is  also 
clear  that  10  samples  per  decorrelation  time  is  suflicient,  at  leasL  for  30  dB  or  less 
fades,  as  long  as  interpolation  is  used  between  realization  points. 

In  summary',  the  adequacy  of  a  realization  depends  on  the  depth  of  fades  that 
are  of  interest.  Realizations  that  are  100  t0  long  accurately  reproduce  the  duration 
and  separation  of  fades  whim  are  less  than  about  20  dB.  If  deeper  fades  are  of  interest, 
then  the  realizations  should  be  longer  and  interpolation  should  be  used. 
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Figure  3.5.  Cumulative  distribution  of  400  r0  realizations  with 
interpolation. 
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Figure  3.0.  Mean  fade  duration  of  400  r0  realizations  with  interpolation. 


ME  AS  FADE  SEPARATION' 


Figure  3.7.  Mean  fade  separation  of  400  r0  realizations  with  interpolation 
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APPENDIX  A 


JOINT  PROBABILITY  DENSITY  FUNCTION  /(«,«') 


The  purpose  of  this  appendix  is  to  derive  the  joint  probability  density  func¬ 
tion  of  the  Rayleigh  amplitude  a  and  its  time  derivative  a1  —  dajdt.  This  function  is 
required  to  calculate  the  temporal  statistics  of  Rayleigh  fading.  A  slightly  different 
form  of  this  derivation  was  first  published  by  Rice  (1948). 

The  starting  point  for  this  calculation  is  the  determination  of  the  joint  prob¬ 
ability  density  function  of  the  in-phasc  and  quadrature  components  x  and  y  of  the 
complex  envelope  of  the  electric  field.  It  is  assumed  that  x  and  y  arc  independent, 
and,  by  the  central  limit  theorem,  that  they  arc  normally  distributed.  Thus  the  joint 
probability  density  function  of  x  and  y  is 


/(*>y) 


1 

2 Vo7 


exp 


x2  4  y2^ 


2o7  \ 


( A  - 1 ) 


Now,  th<'  probability  joint  density  function  of  the  time  derivatives  x'  ~  dx/dt  and 
y'  --  dy/dl  must  be  calculated.  It  will  be  shown  that  x  and  x’  arc  independent,  as 
are  y  and  y‘ .  It  will  be  assumed  that  x,  x1 ,  y,  and  y'  are  jointly  independent.  Thus 
the  joint  probability  density  function  of  x'  and  y'  is  all  that  is  needed  in  addition  to 
Equation  A.l  to  write  down  the  joint  probability  den:  ity  function  /(x,  x',  y,  y') .  Once 
this  function  has  been  obtained,  a  simple  change  of  variables  from  x,x',y,  and  y'  to  a 
and  a  will  yield  the  desired  function. 

In  order  to  determine  the  distribution  of  x'  (or  y'),  consider  the  random 
function  x{l)  written  as  a  Fourier  stochastic  integral 


x(t)  /  z(-jj)  exp(uut)  —  .  (A. 2) 

J  o o  2  71 

The  quantity  2(xi)  is  a  random  function  in  the  Doppler  domain.  It  is  useful  to  assume 
that  x('.u)  is  a  zero-mean,  normally  distributed  random  process,  although  this  is  not 
necessary  because  the  central  limit  theorem  will  make  x(t)  normally  distributed  for 
almost  any  reasonable  distribution  of  a('-a).  However,  with  the  normal  assumption  for 
-j(x'),  x(l)  is  the  sum  of  many  independent,  normally  distributed  random  variables, 
and  is  necessarily  a  zero-menu,  normally  distributed  random  variable. 
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Before  continuing,  it  is  interesting  to  show  the  relationship  between  the  ran¬ 
dom  spectral  components  z(w)  and  the  Doppler  spectrum  ■S'(w).  The  autocorrelation 
function  of  the  stationary  process  x[t)  may  be  written  as 


p[h  ~  ii) 


{ x{ti)x[t2 )) 


d^i  {a(cai)a’(a)2)) 
2rr  a 2 


exp(!WiD  _ 


(A  .3) 


However,  the  correlation  function  p{r)  may  also  be  written  in  terms  of  S (u;) : 


/■» 

S(.o)  exp(fwr)  —  .  (A-4) 

-oo  Z7i 

Note  that,  in  general,  o(-u)  must  be  an  even  function  if  the  autocorrelation  function 
p(r)  is  real.  In  order  to  ensure  that  the  integral  in  Equation  A. 3  is  only  a  function  of 
time  (inference  i  =  t\  -  /2,  the  integrand  must  contain  a  factor  2rr <5 —  ca2).  Using  the 
Dirac  delta  function  to  collapse  the  double  integral  in  Equation  A. 3  and  comparing 
the  result  with  Equation  A. 4  gives 


(.:(w'i)2  (atj),)  'Ino2  —  0^2 ) ■' !>  (at  1)  .  (A. a) 

Tiiis  equation  also  demonstrates  that  the  random  Doppler  spectral  components  of 
z(x)  are  uncorrelaled,  which  is  a  consequence  of  the  assumption  that  the  random 
process  x(t)  is  stationary. 

The  time  derivative  of  i(/j  is  given  by  differentiating  Equation  A. 2,  with  a 
similar  expression  holding  fern  y1: 

(A.(i) 

Because  z ( w )  is  ikh  naliy  distributed  with  zero  mean,  x'{t)  will  also  be  normally 
distributed  with  zero  mean.  The  variance  of  x'[l)  is 


x’{t)  -  /  (rv)a(w)  exi.'(iw/) 

J  CO 
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27 r 
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zix,2\zlxl ).:  (u.’2)/ exppfu-q  -  w2)/i 


<x'(/)/(/) 


A-2 


(A. 7) 
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In  general,  the  variance  of  x'(l)  may  he  written  as 


(x:(<)x'(0) 


( A-8) 


where 


A 


1.0  (Gaussian  Doppler  Spectrum) 

cx/v'l  —  1.5176  (/"‘‘Doppler  Spectrum) 


(A.9) 


and  where  the  parameter  a  is  defined  in  Equation  2.8.  The  cross  correlation  of  x'(t) 
and  x(t)  is 


(x'(i)x(t))  -  f  ^-(t'Ui)  /  ^(z(wl)2*(w2))cxp[i(wi  -  Wj)t! 

J  ■  oo  Lt\  J -ro  in 


dvi 
\n 

-30  du 


io 


r°°  dw 

/  ; — U  6(cu)  =0 

•I  -  «x>  27T 


(A. 10) 


liquations  A. 2  and  A. 5  and  the  fact  that  5(uj)  is  an  even  function  have  been  used 
in  reducing  Equations  A. 7  and  A. 10.  Because  x(t)  and  x'(t)  arc  uncorrclated  and 
normally  distributed,  they  arc  also  independent.  Identical  results  hold  for  the  variance 
of  y1  and  the  cross  correlation  of  y  and  y' . 

The  joint  probability  density  function  of  x,x',y,  and  y'  may  now  be  written 

down: 


A*-*’-**')  -  (^)cxp 


2  2 
2  a2 


WA’)°XP 


*S 


4o2  A2 


.  (A. 11) 


This  function  may  be  transform  to  the  desired  function  of  a  arid  a'  by  making  the 
change  of  variables 


A-3 


x  =  acosO 
y  --  a  sin  6 

The  time  derivatives  of  x  and  y  are 

x1  =  a1  cos  6  -  a6'  sin  6 
y'  —  a  sin  6  i  ad'  cos  6 

which  gives  the  usual  polar  coordinate  equations 

x 2  t  y 2  =  a2 

x'2  +  y'2  =  a'2  +  aV2  . 

The  probability  density  function  coordinate  transformation  is 


(A. 12) 


( A .  13) 


(A.M) 


f(x,  x1  ,y,y')dxdx'dydy'  =  f(a,  a  ,  0,  O')  |  det  ( J)  \dada1  dOdO'  (A.  15) 

where  the  determinate  of  the  Jacobian  of  the  transformation  is 
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(A. 16) 


O'  cos  0 
sin  0 

a'  cos  0  -  aC  sin  0 
a  cos  0 
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A-d 


The  joint  probability  density  function  of  a  and  a'  is  obtained  by  integrating  over  0 
and  O': 


yoo 

f(a,a)  =  /  do  dO'  f  (a,  a' ,  0,  O') 

Jo  J  -oc, 

with  the  result 


(A. 17) 


TO 


27t'/2A  a 


exp 


2  #2 
Tpa 

4A2o2 


(A. 18) 


Thus  it  is  apparent  from  this  equation  that  the  probability  density  function  of  a  is 
Rayleigh;  that  the  probability  density  function  of  a'  is  normal  with  zero  mean  and 
variance  of  2A2o2/rg ;  and  that  a  and  a'  are  independent  because  their  joint  probability 
density  function  is  separable  into  a  function  of  a  times  a  function  of  a' . 


APPENDIX  B 


REALIZATIONS  OF  RAYLEIGH  FADING 


The  methods  of  generating  random  realizations  of  the  impulse  response  func¬ 
tion  for  Rayleigh  fading  channels  are  discussed  in  detail  elsewhere  (Knepp  1982  or 
Dana  198G)  and  arc  only  briefly  reviewed  in  this  appendix.  The  starting  point  of 
the  method  of  generating  a  realization  of  flat  Rayleigh  fading  is  the  Doppler  Power 
Spectral  Density  (PSD)  function,  5(ca).  This  function  has  two  commonly  used  forms 
which  were  discussed  in  Section  2: 


and 


S(u>)  -  \/7r r0  exp 


T2'  .2 
T0'° 


(Gaussian  PSD) 


(IU) 


4r0  1 

Q  |l  i  rjca2/a2j2 


(D4  PSD) 


(P.2) 


where  a  -  2.146193  in  the  f " 4  form  of  the  PSD.  The  quantity  f?(ta)dcu/'27r  is  the  mean 
fraction  of  power  in  the  Doppler  frequency  interval  u>/ 2n  to  (u  t  (Lj)/2n. 

The  discrete  realizations  of  the  channel  impulse  response  function  will  con¬ 
tain  A’  time  samples  with  A’0  samples  per  decorrelation  time.  Thus  the  time  spacing 
of  tin;  discrete  samples  is 


At  = 


Tq 


(11.3) 


am!  the  total  time  duration  of  the  realization  is  NAt.  In  the  Doppler  frequency 
domain,  the  spacing  of  the  discrete  samples  is 


A-j 


2  71 

X~A t 


(ll,t) 


Note  that  the  quantity  A wAl,  which  will  appear  later  in  a  Fourier  trans  arm,  is  equal 
to  2 a  /  A; . 


P-1 


The  samples  in  the  frequency  domain  are  generated  by  first  calculating  the 
fraction  of  signal  power  in  each  Doppler  frequency  bin,  S}  —  S(jAuj)  Aw/27r.  For  the 
Gaussian  PSD, 


jvv|i 

N‘  j 

and  for  the  /'  4  PSD, 


s,  --=  -v 


V*A'o 

N~ 


exp 


(J  -  -N/2,  •  •  • , yV/2  --  1)  (H.5) 


_  4A'0 

lj  aN 


jW/V2 

1  |  - - - 

a2/V2 


[j  =  --N/2,  -  ,  N/2  —  1)  (15.(3) 


Next,  the  random  Doppler  frequency  spectrum  //(jAw)  of  the  impulse  response  func¬ 
tion  is  generated: 


11  (jAw)  =  2^\/%  •  (»-7) 

The  leading  factor  27t/Acj  has  been  included  so  that  the  discrete  Fourier  transform 
of  II(jA'jj)  will  be  dimensionless.  'File  random  components  of  the  spectrum,  £;,  are 
complex,  normally  disturbed  random  variables  with  the  properties: 


(66)  =-  <5;.*  (»•») 

(66)  =  b  , 

where  the  quantity  6a  is  the  Kronecker  delta  symbol,  'fhe  random  samples  of  £  may 
be  generated  using 


6  -  \/ "M'Mj)  exp(27uu2i;)  (13.9) 

where  Uj  ;  and  u2 ;  are  independent  random  variables  uniformly  distributed  on  the 
interval  (0, 1 ) . 


Finally,  the  random  Doppler  spectrum  of  the  channel  impulse  re. mouse  func¬ 
tion  is  Fourier  transformed  to  the  time  domain.  In  continuous  notation,  us  Fourier 
transform  is 
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r  co  < 

h(t )  -  /  11  (  ''exp(twi); 

J  —  'VJ 


and  in  discrete  notation 


"/2-1  Aw 

h(kAt)  =  // (j  Aw)  exp[r(j  Aw)  (/: At)  j  — - 

i—~NJ2 

N/2-l 

=  \fs~)Zjcxp\t2nJk/N\  ' 
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where  A:  =  0,  1 ,  •  •  • ,  iV  —  1. 
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